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Abstract 

Background: The associations between Rad51 gene polymorphisms (G135C and G172T) and risk of cancer have been 
investigated, but the results were inconclusive. To get a comprehensive evaluation of the association above, we performed 
a meta-analysis of published studies. 

Methods: A computerized search of PubMed, Embase and Web of Knowledge databases for all relevant studies was 
performed and the data were analyzed in a meta-analysis. The overall odds ratio (OR) with the 95% confidence interval (95% 
CI) was calculated to assess the strength of the association between Rad51 polymorphisms and cancer risk. Data were 
analyzed using fixed- or random-effects model when appropriate. Sensitivity analysis and publication bias test were also 
estimated. 

Results: Overall, a total of 54 case-control studies were included in the current meta-analysis, among which 42 studies with 
19,142 cases and 20,363 controls for RAD51 G135C polymorphism and 12 studies with 6,646 cases and 6,783 controls for 
G172T polymorphism. For G135C polymorphism, the pooled results indicated that significantly increased risk was found in 
overall cancers (homozygote model: OR = 1.776, 95% CI = 1 .288-2.449; allelic genetic model: OR =1.169, 95% CI = 1.016- 
1.345; recessive model: OR = 1.946, 95% CI = 1.336-2.835), especially in breast cancer (homozygote model: OR = 1.498, 95% 
CI = 1 .026-2.1 89; recessive model: OR = 1 .732, 95% CI = 1 .1 70-2.562). For G1 72T polymorphism, a decreased cancer risk was 
observed in head and neck cancer (homozygote model: OR = 0.621, 95% CI = 0.460-0.837; allelic genetic model: OR = 
0.824, 95% CI = 0.716-0.948; recessive model: OR = 0.639, 95% CI = 0.488-0.837). 

Conclusions: Our results suggested that the Rad51 G135C polymorphism is a candidate for susceptibility to overall cancers, 
especially to breast cancer, and that the Rad51 G172T might play a protective role in the development of head and neck 
cancer. 
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Introduction 

Human cancer is still one of the leading causes of death 
worldwide, resulting in one of the most challenging global health 
issues confronted by mankind today. According to etiological 
studies, carcinogenesis of cancer is a complex, multistep and 
multifactor process, in which many genetic and environmental 
factors are involved. In recent years, it has become clear that 
individual variation in genetic backgrounds can lead to various 
consequences following the environmental exposure and may 
ultimately contribute to the cancer pathogenesis and progression 
[1-3]. 

DNA repair pathways are responsible for maintaining the 
genomic stability and integrity and play a pivotal role in protecting 
against genetic mutations [4]. DNA repair genes have been 



proposed as considerable factors in the prevention of genomic 
damage and continuously monitor chromosomes to correct 
injuries caused by exogenous agents such as ultraviolet light or 
cigarette smoke, and endogenous mutagens [5,6]. Recent reports 
have indicated that genetic variation in DNA repair genes could 
cause altered DNA repair capacity, leading to accumulation of 
DNA damage, followed by programmed cell death or unregulated 
cell growth and may account, in part, for the cancer development 
[7]- 

Human RAD51, one of the key proteins for homologous 
recombination, is essential to meiotic and mitotic recombination 
and plays a crucial role in homologous recombination repair of 
DNA double-strand breaks [8]. It functions by forming nucleo- 
protein filaments on single stranded DNA, inducing homologous 
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pairing and mediating strand exchange reactions between single 
and double stranded DNA during repair [9] . The RAD5 1 gene is 
located on chromosome 15ql5.1 in humans and thought to 
participate in a common double-strand break repair pathway. In 
recent years, RAD51 gene polymorphisms have attracted wide- 
spread attention. Two commonly studied polymorphisms of 
RAD51 gene are G135C (rs 180 1320), a G to C transversion at 
position +135, and G172T (rsl801321), a G to T transversion in 
the 172 position. Both of them are located in the 5' untranslated 
region and seem to be of functional relevance. These two 
polymorphisms were shown to affect mRNA stability or transla- 
tional efficiency, leading to altered polypeptide product levels and 
altering the function of encoding RAD5 1 protein, and influenced 
the DNA repair capacity to some extent [10,11]. 

In the past decade, a number of molecular epidemiological 
studies have been done to evaluate the association between 
RAD51 gene polymorphisms (G135C and G172T) and cancer risk 
in diverse populations, but the results remained controversial. 
Therefore, to derive a more precise estimation of the association 
between RAD51 G135C and G172T polymorphisms and cancer 
risk, a meta-analysis was performed. To the best of our knowledge, 
this is the most comprehensive meta-analysis regarding the 
RAD51 polymorphisms and cancer risk. 

Materials and Methods 

Search strategy and data extraction 

All studies investigating the association between the RAD51 
gene polymorphisms (G135C and G172T) and risk of cancer were 
identified by comprehensive computer-based searches of PubMed, 
Embase and Web of Knowledge databases (the last search update 
on August 25, 2013). The search was performed using various 
combinations of keywords like ("RAD51 gene" OR "RAD51 
recombinase gene") AND ("polymorphism" OR "variant" OR 
"variants"). The exact search is available on request from the 
authors. Additional studies were also identified by a hand search of 
all the references of retrieved articles. Our search was restricted to 
studies published in the English language. 

Inclusion and Exclusion criteria 

Studies included in the current meta-analysis had to meet all the 
following criteria: (1) studies to investigate the associations between 
the polymorphisms of G135C or G172T in RADSlgene and risk 
of cancer; (2) an unrelated case-control or cohort design (3) 
sufficient data (genotype distributions for cases and controls) to 
calculate an odds ratio (OR) with its 95%CI; (4) studies published 
in English; (5) genotype distribution of control population 
consistent with Hardy-Weinberg Equilibrium (HWE). We did 
not consider abstracts or unpublished reports. Case reports, 
editorials, review articles, and letters were excluded. Articles were 
also excluded if they did not include a control population and did 
not determine genotype frequency. If studies with the same or 
overlapping data were published by the same authors, the study 
with the larger sample size was selected. The supporting PRISMA 
checklist is available as supporting information; see Checklist S 1 . 

Data extraction 

Two authors extracted information from all eligible publications 
independendy according to the inclusion criteria listed above. 
Disagreement was resolved by the evaluation of a third reviewer 
and discussion until a consensus was reached. The following 
characteristics were collected from each study: the first author, 
year of publication, country, patient ethnicity, cancer type, source 
of control groups (population- or hospital-based controls or mixed 



(composed of both population- and hospital-based controls)), and 
genotype frequencies in case and control groups. Meanwhile, we 
did not define any minimum number of cases or controls to be 
included in our meta-analysis. 

Statistical analysis 

We first analyzed HWE in the controls for each study using 
goodness-of-fit test (chi-square or Fisher's exact test) and violation 
of HWE was determined by P<0.05. Crude odds ratios (ORs) 
with 95% confidence intervals (CIs) were used to assess the 
strength of association between the RAD5 1 gene polymorphisms 
and cancer susceptibility. The pooled ORs for RAD51 G135C 
polymorphism were performed under dominant model (CC+GC 
vs. GG), recessive model (CC vs. GG+GC), homozygote model 
(CC vs. GG) and allelic genetic model (C vs. G). C and G 
represent the minor and the major allele respectively. The same 
methods were applied to the analysis of the RAD51 G172T 
polymorphism. Stratified analyses were conducted with respect to 
ethnicity, cancer type and source of controls. 

A % 2 -based CHest was performed to test the heterogeneity 
across the eligible comparisons, which is considered to be 
significant if P<0.05. The variation caused by heterogeneity was 
estimated by calculating the inconsistency index I 2 , with I 2 <25%, 
25-75% and >75% representing low, moderate or high degrees of 
inconsistency, respectively [12]. The pooled OR was calculated by 
a fixed-effects model (the Mantel-Haenszel method) if the result of 
the 0_ test was P>0.05, which indicated that the between-study 
heterogeneity was not significant [13]. Otherwise, a random- 
effects model (the Der-Simonian and Laird method) was used [14]. 
Sensitivity analysis was carried out by removing each study at a 
time to evaluate the stability of the results under either genotypic 
models or the allelic model. Additionally, Begg's test and Egger's 
linear regression test by visual inspection of the funnel plot were 
carried out to address the potential publication bias and P<0.05 
was considered as an indicator of significant publication bias [15]. 

All statistical analyses were performed using the STATA 
software (version 11; Stata Corporation, College Station, Texas). 
Two-sided P values less than 0.05 were considered significant. 

Results 

Studies included in the meta-analysis 

The initial literature search through PubMed, Embase and Web 
of Knowledge databases yielded 203 published articles after 
duplicates were removed. When reviewed by tides or abstracts, 
1 1 5 records did not fulfill the inclusion criteria, leaving 88 
potentially relevant studies that were reviewed in full-text. Among 
the remaining 88 articles, 10 were not concerned with G135C or 
G172T polymorphisms in RAD51gene, 7 were not human studies, 
4 was not published in English, 6 were not case-control studies, 5 
were no usable reported data, 2 were meeting abstracts, 4 were 
meta-analysis, and 1 1 were not in HWE; these publications were 
also excluded. Finally, a total of 54 case-control studies in 37 
articles were identified in the current meta-analysis [16-54], 
among which 42 with 19142 cases and 20363 controls for RAD51 
G135C polymorphism and 12 with 6646 cases and 6783 controls 
for G172T polymorphism. Genotype distributions in the controls 
of all selected studies are in agreement with HWE. The flow of 
study selection is shown in Figure 1, and the main characteristics 
of eligible studies were shown in Table 1 and Table 2. 
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203 articles and related reference indentified 
through database searching 





115 Excluded 




43 were obviosly irrelevant 




21 not for cancer research 




27 not on Rad51 gene 




17 not on human studies 




7 reviews 



88 abstracts considered for further evaluation 



31 Excluded 
10 not related to these two polymorphisms 
7 not human studies 
4 meta-analysis 
4 not published in English 
6 not case-control studies 



57 of full-text articles assessed for eligibility 



18 Excluded 
11 not in HWE 
5 no sufficient usable date 
2 meeting abstracts 



39 articles included in analysis 



Figure 1. The flow diagram of the literature search and the study selection. 

doi:10.1371/journal.pone.0087259.g001 



Meta-analysis result 

The pooled results of meta-analysis for the association between 
RAD51 polymorphisms (G135C and G172T) and cancer suscep- 
tibility are shown in Tables 3 and Tables 4. 

As for G135C polymorphism, a total of 42 case-control studies 
in 37 publications with 19,142 cases and 20,363 controls were 
identified. Overall, significandy elevated cancer risk was found in 
all genetic models (homozygote model: OR= 1.776, 95% 
01= 1.288-2.449, Figure 2; allelic genetic model: OR =1.169, 
95% CI = 1.016-1.345; recessive model: OR= 1.946, 95% 
CI= 1.336-2.835) except in dominant model (OR= 1.039, 95% 
CI = 0.942-1.146). The heterogeneity was significant in all genetic 
models and the detailed data are shown in Table 3. These eligible 
studies were analyzed by stratified analysis. In the stratified 
analysis of the effect of cancer types, a significant association was 
found for breast cancer (homozygote model: OR= 1.498, 95% 
CI= 1.026 -2.189; recessive model: OR = 1.732, 95% 
CI = 1.170-2.562). However, no significant association with 
cancer risk was demonstrated in overall population with ovarian 
cancer, colorectal cancer, acute myelocytic leukemia as well as 
head and neck cancers. As for ethnicity, our results showed G135C 
polymorphism was associated with increased risk of cancer among 
all populations under homozygote model and recessive model. 
When stratified based on source of controls, significantly increased 
risks were also observed in both population-based control 
subgroups and hospital-based control subgroups (Table 3). 



With respect to G172T polymorphism, a total of 12 case-control 
studies in 9 publications with 6,646 cases and 6,783 controls were 
selected. As shown in Table 4, the pooled results revealed no 
significant associations between G172T polymorphism and cancer 
susceptibility in all genetic models (homozygote model: 
OR= 1.014, 95% CI = 0.872-1. 173; dominant model: 
OR= 0.980, 95% CI = 0.906-1.061, Figure 3; recessive model: 
OR= 1.011, 95% CI = 1.241-14.879; allelic genetic model: 
OR= 0.993, 95% CI = 0.941-1.048). The heterogeneity was 
significant in all genetic models except for dominant model 
(P = 0.414). We also analyzed these eligible studies by stratified 
analysis. As we divided the studies by cancer type, the result 
suggested that a decreased cancer risk was found in head and neck 
cancers (homozygote model: OR = 0.621, 95% CI = 0.460-0.837; 
allelic genetic model: OR = 0.824, 95% CI = 0.716-0.948; reces- 
sive genetic model: OR = 0.639, 95% CI = 0.488-0.837) Never- 
theless, we did not find significant association between G172T 
polymorphism and breast cancer and ovarian cancer. When 
stratified according to ethnicity, the result showed no evidence that 
the G172T polymorphism was significantly associated with an 
increased cancer risk in Caucasian populations. In the subgroup 
analysis by source of controls, no significant association with 
cancer risk was observed in both population-based and hospital- 
based control subgroups (Table 4). 
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Table 1. Characteristics of the studies included on Rad51 G135C polymorphism. 





First author 


Year 


Country 


Ethnicity 


Cancer type 


Source of 
control 


Cases 


Controls 


GG 


GC 


CC 


GG 


GC 


CC 


HWE 


Kuschel 


2002 


Germany 


Caucasian 


Breast 


PB 


2172 


840 


1904 


255 


13 


722 


116 


2 


0.23 


Seedhouse 


2003 


UK 


Caucasian 


AML 


PB 


186 


206 


166 


18 


2 


171 


32 


3 


0.30 


Wang 


2004 


USA 


Caucasian 


Glioma 


Mixed 


309 


342 


265 


40 


4 


301 


41 


0 


0.24 


Webb (1) 


2005 


Australia 


Australian, Caucasian 


Breast 


PB 


1444 


788 


1221 


212 


11 


676 


104 


8 


0.08 


Webb (2) 


2005 


Australia 


Australian, Caucasian 


Ovarian 


PB 


546 


1126 


457 


85 


4 


971 


145 


10 


0.08 


Auranen (1) 


2005 


UK 


Caucasian 


Ovarian 


PB 


729 


847 


642 


84 


3 


745 


100 


2 


0.48 


Auranen (2) 


2005 


UK 


African 


Ovarian 


PB 


326 


419 


270 


52 


4 


357 


61 


1 


0.24 


Auranen (3) 


2005 


UK 


Caucasian 


Ovarian 


PB 


278 


699 


241 


36 


1 


616 


78 


5 


0.15 


Auranen (4) 


2005 


UK 


Caucasian 


Ovarian 


PB 


296 


840 


266 


29 


1 


722 


116 


2 


0.23 


Lee 


2005 


Korea 


Asian 


Breast 


HB 


782 


587 


611 


143 


28 


450 


123 


14 


0.11 


Sliwinski 


2005 


Poland 


Caucasian 


Breast 


HB 


150 


150 


108 


38 


4 


106 


41 


3 


0.67 


Hanna 


2005 


Poland 


Caucasian 


Breast 


PB 


100 


106 


31 


40 


29 


21 


48 


37 


0.45 


Romanowicz 


2006 


Poland 


Caucasian 


Breast 


PB 


100 


106 


31 


40 


29 


21 


48 


37 


0.45 


Tarasov 


2006 


Russia 


Caucasian 


Breast 


PB 


151 


191 


111 


36 


4 


148 


41 


2 


0.65 


Rollinson 


2006 


UK 


Caucasian 


AML 


HB 


466 


936 


431 


34 


1 


817 


115 


4 


0.98 


Antoniou 


2007 


UK 


Caucasian, African 


Breast 


HB 


4443 


4069 


3838 


567 


38 


3485 


565 


19 


0.44 


Costa 


2007 


Portugal 


Caucasian 


Breast 


PB 


265 


435 


216 


45 


4 


381 


53 


1 


0.55 


Lu 


2006 


USA 


Caucasian 


HNC 


HB 


716 


719 


624 


91 


1 


622 


96 


1 


0.17 


Jakubowska (1) 


2007 


Poland 


Caucasian 


Breast 


HB 


258 


258 


210 


48 


0 


188 


68 


2 


0.12 


Jakubowska (2) 


2007 


Poland 


Caucasian 


Ovarian 


HB 


127 


127 


104 


23 


0 


89 


37 


1 


0.17 


D.Figueroa 


2007 


USA 


Caucasian 


Bladder 


HB 


1085 


1032 


932 


147 


6 


909 


116 


7 


0.12 


Voso 


2007 


Italy 


Caucasian 


AML 


HB 


160 


161 


125 


33 


2 


142 


18 


1 


0.61 


Jara 


2007 


Chile 


Latin-American 


Breast 


PB 


131 


247 


113 


16 


2 


222 


25 


0 


0.40 


Werbrouck 


2008 


Belgium 


Caucasian 


HNC 


HB 


152 


157 


136 


15 


1 


134 


23 


0 


0.32 


Synowiec 


2008 


Poland 


Caucasian 


Breast 


PB 


41 


48 


18 


10 


13 


17 


27 


4 


0.14 


Jakubowska 


2008 


Poland 


Caucasian 


Breast 


PB 


1007 


1069 


785 


207 


15 


822 


232 


15 


0.76 


Bhatla 


2008 


USA 


Caucasian, African 


AML 


HB 


452 


646 


374 


73 


5 


555 


85 


6 


0.18 


Krupa 


2009 


Poland 


Caucasian 


Breast 


PB 


135 


175 


91 


33 


11 


105 


63 


7 


0.52 


Dhillon 


2009 


Australia 


Caucasian 


Prostate 


HB 


116 


132 


97 


18 


1 


119 


13 


0 


0.55 


Jara 


2010 


Chile 


Latin-American 


Breast 


PB 


267 


500 


232 


33 


2 


441 


58 


1 


0.53 


Krupa 


2010 


Poland 


Caucasian 


Endometrial 


PB 


30 


30 


6 


8 


16 


19 


9 


2 


0.52 


Romanowicz 


2010 


Poland 


Caucasian 


Breast 


HB 


220 


220 


141 


69 


10 


157 


58 


5 


0.89 


Liu 


2011 


China 


Asian 


AML 


HB 


105 


704 


72 


25 


8 


511 


175 


18 


0.52 


Gil 


2011 


Poland 


Caucasian 


Colorectal 


HB 


133 


100 


100 


29 


4 


73 


27 


0 


0.12 


Hamdy 


2011 


Egypt 


Caucasian 


AML 


PB 


50 


30 


39 


9 


2 


26 


3 


1 


0.06 


Romanowicz 


2012 


Poland 


Caucasian 


Endometrial 


HB 


230 


236 


40 


25 


165 


59 


132 


45 


0.06 


Gresner 


2012 


Poland 


Caucasian 


HNC 


PB 


81 


87 


67 


13 


1 


71 


14 


2 


0.22 


Zhang 


2012 


China 


Asian 


Cervical 


HB 


80 


175 


58 


20 


2 


122 


50 


3 


0.41 


Mucha 


2012 


Poland 


Caucasian 


Colorectal 


HB 


200 


200 


161 


34 


5 


157 


37 


6 


0.05 


Romanowicz 


2012 


Poland 


Caucasian 


Colorectal 


HB 


320 


320 


51 


56 


213 


91 


164 


65 


0.57 


Romanowicz 


2012 


Poland 


Caucasian 


HNC 


PB 


253 


253 


174 


69 


10 


190 


58 


5 


0.82 


Smolarz 


2013 


Poland 


Caucasian 


Breast 


PB 


50 


50 


8 


8 


34 


14 


26 


10 


0.74 



PB: population based; HB: hospital based; HWE: Hardy-Weinberg equilibrium (significant at the 0.05 level} 
AML: acute myelocytic leukemia; HNC: head and neck cancer. 
doi:1 0.1 371 /journal.pone.0087259.t001 



Sensitive analysis 

Given that the significant between-study heterogeneity for 
RAD51 G135C and G172T polymorphisms, the random-effect 
model was used to calculate the pooled results if the heterogeneity 



was significant. Meanwhile, we also performed a sensitivity 
analysis to assess the effects of each study on the pooled ORs by 
omission of individual studies. The sensitivity analysis showed that, 
for each polymorphism, no single study qualitatively changed the 
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Table 2. Characteristics of the studies included on Rad51 G172T polymorphism. 





First author 


Year 


Country 


Ethnicity 


Cancer 
type 


Source of 
control 


Cases 


Controls 


GG 


GC 


TT 


GG 


GC 


TT 


HWE 


Kuschel 


2002 


Germany 


Caucasian 


Breast 


PB 


2235 


736 


744 


1061 


430 


226 


371 


139 


0.54 


Auranen(l) 


2005 


UK 


Caucasian 


Ovarian 


PB 


730 


847 


226 


363 


141 


273 


433 


141 


0.16 


Auranen(2) 


2005 


UK 


African 


Ovarian 


PB 


321 


412 


119 


145 


57 


149 


189 


74 


0.30 


Auranen(3) 


2005 


UK 


Caucasian 


Ovarian 


PB 


293 


607 


112 


130 


51 


235 


277 


95 


0.38 


Auranen(4) 


2005 


UK 


Caucasian 


Ovarian 


PB 


300 


736 


94 


157 


49 


226 


371 


139 


0.54 


Lee 


2005 


Korea 


Asian 


Breast 


HB 


784 


591 


721 


54 


9 


533 


54 


4 


0.05 


Rollinson 


2006 


UK 


Caucasian 


AML 


HB 


469 


940 


144 


225 


100 


331 


445 


164 


0.49 


Lu 


2006 


USA 


American 


HNC 


HB 


716 


719 


261 


351 


104 


240 


335 


144 


0.17 


Silva 


2009 


Portugal 


Caucasian 


Breast 


HB 


288 


548 


94 


139 


55 


168 


275 


105 


0.69 


Gresner 


2012 


Poland 


Caucasian 


HNC 


PB 


81 


110 


36 


43 


2 


43 


54 


13 


0.52 


Romanowicz 


2012 


Poland 


Caucasian 


Colorectal 


HB 


320 


320 


81 


150 


89 


84 


142 


94 


0.05 


Bastos 


2009 


Portugal. 


Caucasian 


Thyroid 


HB 


109 


217 


28 


51 


30 


76 


98 


43 


0.27 



PB: population based; HB: hospital based; HWE: Hardy-Weinberg equilibrium (significant at the 0.05 level} 
AML: acute myelocytic leukemia; HNC: head and neck cancer 
doi:1 0.1 371 /journal.pone.0087259.t002 



pooled ORs, suggesting that the results of this meta-analysis were 
statistically stable and reliable (Figure_Sl and Figure _S2). 

Publication bias diagnostics 

We further identify the potential publication biases of literatures 
by Egger's test and funnel plot. In all studies, no funnel plot 
asymmetry was found. The results of the Egger's test for RAD5 1 
G135C and G172T polymorphisms did not show any evidence of 
publication bias (t = — 1.11, P = 0.275 for G135C under homozy- 
gote comparison model, Figure 4; t = -0.09, P = 0.929 for G 1 72T 
under homozygote comparison model, Figure 5). 

Discussion 

It is well reported that double-strand break damage is the most 
dangerous lesion observed in eukaryotic cells because it may cause 
cell death or constitute a serious threat to cell viability and genome 
stability. It has the potentiality to permanently arrest cell cycle 
progression and endanger cell survival [55]. Since DNA repair 
mechanisms are essential to preserve genomic stability and 
functionality, defects in DNA repair can result in the development 
of chromosomal aberrations which may lead to an increased 
susceptibility to cancer [4,56,57]. Homologous recombination and 
non-homologous end joining have been extensively studied as two 
distinct pathways in the repair of double-strand breaks in 
mammalian cells. Homologous recombination is a high-fidelity 
process that utilizes DNA sequence, a sister chromatid or 
homologous chromosome in close proximity to the break as a 
template [58-60]. In this repair process, an early procedure is the 
resection of the 3'ends of the DSBs to form single stranded tails 
that invade the intact homologous DNA double helix forming a 
Holliday junction [61,62]. RAD51, a kind of ubiquitous strand 
exchange protein, is known to be a central enzyme involved in 
DNA double-strand break repair by homologous recombination. 
It could polymerize onto single-stranded DNA and searches for 
homology in a duplex donor DNA molecule, usually the sister 
chromatid [63]. Recent researches have suggested two common 
polymorphisms (G135C and G172T) located in the 5' untrans- 
lated region seems to be of functional relevance. Furthermore, 
many functional studies revealed that these polymorphisms could 



affect mRNA stability or translational efficiency, resulting in 
changes in both polypeptide product levels and the function of 
encoding RAD51 protein, and thus influenced the DNA repair 
capacity to some extent [10,11]. In addition, the association of 
Rad51 variants (G135C and G172T) and risk of cancer has been 
extensively investigated in different populations. However, the 
results of these studies were inconsistent. Therefore, we conducted 
a meta-analysis to summarize the effects of Rad5 1 variation on risk 
of cancer. 

In this meta-analysis, 54 case-control studies (42 for G135C 
polymorphism, 12 for G172T polymorphism) were performed to 
provide the most comprehensive assessment of the relationship 
between RAD51 polymorphisms and cancer risk. For Rad51 
G135C polymorphism, the C allele of G135C polymorphism had 
significant association with the cancer susceptibility for the 
homozygote model, allelic genetic model, and recessive genetic 
model in overall populations. Nevertheless, the results suggested 
that Rad51 G172T polymorphism was not associated with overall 
cancer risk when all studies were accumulated together. Consid- 
ering the possible role of ethnic differences in genetic backgrounds, 
we performed subgroup analysis based on ethnicity. Consequently, 
significant association was found in both Caucasians and Asians 
for Rad51 G135C polymorphism but not for G172T polymor- 
phism. When stratified by the source of controls, our results found 
evidence of an association between cancer risk and G135C 
polymorphism in both population-based and hospital-based 
controls, while no significant association was indicated in either 
population-based or hospital-based controls for G172T polymor- 
phism. In the stratified analysis by cancer type, our results strongly 
indicated that Rad51 G135C polymorphism was associated with 
increased breast cancer risk while G172T polymorphism with 
decreased head and neck cancer risk. 

Previous meta-analyses were carried out to assess the effect of 
Rad51 G135C polymorphism on either the risk of breast cancer or 
acute leukemia [64,65] . Comparing with them, our study has some 
improvements. First, this is the first report not only to analyze two 
polymorphisms in Rad51 gene (G135C and G172T) and cancer 
risk in different cancer forms, but also to identify the G172T 
polymorphism as a risk factor for head and neck cancers. Second, 
we provided a more comprehensive data analysis by calculating 
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Figure 2. Forest plot for association of Rad51 G135C polymorphism and cancer risk (homozygote model, CC vs. GG). 

doi:1 0.1 371 /journal.pone.0087259.g002 



four different genetic models and performing subgroup analysis 
based on ethnicity, cancer types and source of controls. Third, we 
excluded those studies in which genotype distributions in the 
controls were not in agreement with HWE because they could 
influence the results. 

Heterogeneity between studies should be noted because it may 
affect the strengths of the meta-analysis. In the current meta- 
analysis, significance heterogeneity was observed for both Rad5 1 
G135C and G172T polymorphisms. Thus, random-effect models 
were used if significant heterogeneity was identified. Meanwhile, to 
diminish the heterogeneity, we carried out subgroup analysis based 
on ethnicity, cancer types and source of controls. The results 
indicated that heterogeneity reduced or disappeared in subgroups. 
We also performed sensitivity analysis to ascertain the primary 
origin of the heterogeneity. The analysis showed that no single 



study materially altered the pooled ORs, suggesting that the results 
of this meta-analysis were statistically stable and reliable. The 
publication bias for the association between these two polymor- 
phism and cancer risk was not observed. 

Some limitations of the present meta-analysis should be taken 
into consideration when interpreting the results. First of all, only 
published studies and papers written in English were searched in 
this meta-analysis, some unpublished studies or studies written in 
other language that might also meet the inclusion criteria were 
overlooked. Second, in some studies, detailed information such as 
age and sex in case and control of different genotypes were not 
available, which limited further estimates to a certain extent. 
Third, the current meta-analysis did not consider gene-gene and 
gene-environment interactions due to the lack of sufficient data. 
Further studies are needed to evaluate the possible gene-gene and 
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Figure 3. Forest plot for association of Rad51 G172T polymorphism and cancer risk (dominant model, TT+GT vs. GG). 

doi:1 0.1 371 /journal.pone.0087259.g003 



gene-environment interactions in the association between Rad51 
gene polymorphism and susceptibility to cancer. Fourth, most of 
the patients in our study were Caucasians, which may limit the 
general application of our results. In spite of these, our present 
meta-analysis also had some advantages. First, analyzing two 



Rad5 1 gene polymorphisms with a total of 54 case-control studies 
has a much greater statistical power compared with any single 
study. Second, we excluded the studies in which genotype 
frequencies in controls were not in accordance with HWE, 
providing sufficient evidence for drawing safe conclusions about 



Begg's funnel plot with pseudo 95% confidence limits 
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Figure 4. Begg's funnel plot for publication bias in studies on Rad51 G135C polymorphism and cancer (homozygote model, CC vs. 
GG). 

doi:10.1371/journal.pone.0087259.g004 
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Begg's funnel plot with pseudo 95% confidence limits 
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Figure 5. Begg's funnel plot for publication bias in studies on Rad51 G172T polymorphism and cancer (homozygote model, TT vs. 
GG). 

doi:1 0.1 371 /journal.pone.0087259.g005 



the association between the Rad51 polymorphisms and cancer 
risk. Third, the stability and credibility of the present meta-analysis 
was confirmed by the sensitivity analyses and publication biases 
analyses. Last, the findings highlight the association between 
Rad51 gene polymorphisms and cancer development and will 
provide directions for future research on molecular mechanism of 
cancer. 

Conclusions 

Our investigations suggested that the Rad51 G135C polymor- 
phism is a candidate for susceptibility to overall cancers, especially 
to breast cancer, and that the G172T polymorphism is signifi- 
cantly associated with decreased risk of head and neck cancers. 
Further studies are needed with large sample size and deeper 
evaluation about the effect of gene-gene and gene-environment 
interactions on the Rad5 1 polymorphisms and cancer risk. 
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